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Abstract 

The first phase of stellar evolution in the history of the Universe 
may be Dark Stars, powered by dark matter heating rather than by 
nuclear fusion. Weakly Interacting Massive Particles, which may be 
their own antipartners, collect inside the first stars and annihilate to 
produce a heat source that can power the stars for millions to billions 
of years. In this paper we show that these objects can grow to be 
supermassive dark stars (SMDS) with masses ^ (10^ — 10'^)Mq. The 
growth continues as long as dark matter heating persists, since dark 
stars are large and cool (surface temperature ^ 5 x lO^K) and do not 
emit enough ionizing photons to prevent further accretion of baryons 
onto the star. The dark matter may be provided by two mechanisms: 
(1) gravitational attraction of dark matter particles on a variety of 
orbits not previously considered, and (2) capture of WIMPs due to 
elastic scattering. Once the dark matter fuel is exhausted, the SMDS 
becomes a heavy main sequence star; these stars eventually collapse 
to form massive black holes that may provide seeds for supermassive 
black holes in the Universe. SMDS are very bright, with luminosities 
exceeding (10^ — W^^)Lq. We demonstrate that for several reasonable 
parameters, these objects will be detectable with JWST. Such an ob- 
servational discovery would confirm the existence of a new phase of 
stellar evolution powered by dark matter. 
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1 Introduction 



Spolyar et al. (2008) first considered the effect of dark matter (DM) particles 
powering the first stars. The first stars formed when the universe was about 
200 milhon years old, at z = 10 — 50, in IO^Mq haloes consisting of 85% DM 
and 15% baryons in the form of H and He from big bang nucleosynthesis; 
for reviews of the standard picture of the formation of the first stars, see 
Barkana & Loeb (2001), Yoshida et al. (2003), Bromm & Larson (2004), 
and Ripamonti & Abel (2005). The canonical example of particle DM is 
Weakly Interacting Massive Particles (WIMPs). In many theories WIMPs 
are their own antiparticles and annihilate among themselves wherever the 
DM density is high. Recently there has been much excitement in the dark 
matter community about possible detections of WIMPs via annihilation to 
positrons seen by the PAMELA satellite (Adriani et al. 2009); annihilation 
to 7-rays seen by FERMI (Abdo et al. 2009, Dobler et al. 2009), and in 
the direct detection experiments DAMA and CDMS (Bernabei et al. 2009, 
Ahmed et al. 2009). The annihilation rate is n'^{av) where is WIMP 
density and we take the standard annihilation cross sectior0 

{av) = 3 X lO^^^cmVs, (1) 

and WIMP masses in the range 1 GeV-10 TeV. The first stars are particularly 
good sites for annihilation: they form in the right place — in the high density 
centers of DM haloes — and at the right time — at high redshifts (density 
scales as (1 + z)^). Spolyar et al. (2008) found that above a certain density 
(«i 10^^ cm~^) the WIMP annihilation products remain trapped in the star, 
thermalize with the star, and thereby provide a heat source. These first 
Dark Stars (DS) are stars made primarily of hydrogen and helium with only 
~ 0.1% of the mass in the form of DM; yet they shine due to DM heating. 
Note that the term 'Dark' refers to the power source, not the appearance or 
the primary matter constituent of the star. 

Dark stars are born with masses ~ IMq. They are giant puffy (~ 10 
AU), cool (surface temperatures < 10,000K), yet bright > 10^ Lq objects 
(Freese et al. 2008a). They reside in a large reservoir (~ 10^ Mq) of baryons, 
i.e., ~ 15% of the total halo mass. These baryons can start to accrete onto 
the dark stars. Previous work (Freese et al. 2008a; Spolyar et al. 2009) 

^Annihilation in the early universe with this value of the cross section leaves behind the 
correct relic WIMP DM abundance today, ^ 24% of the energy density of the universe. 
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followed the evolution of dark stars from their inception at IMq, as they 
accreted baryons from the surrounding halo, up to ~ lOOOM©. Dark stars 
can continue to grow in mass as long as there is a supply of DM fuel. 

The exciting new development of this paper is that we follow the growth 
of dark stars to become supermassive dark stars (SMDS) of mass M* > 
IO^Mq. Specifically we study the formation of lO^M© SMDS in lO^M© DM 
haloes and 10'' SMDS in IO^Mq haloes; perhaps the SMDSs become even 
larger. Hoyle & Fowler (1963) first postulated the existence of such large stars 
but were not aware of a mechanism for making them. Now the confluence 
of particle physics with astrophysics may be providing the answer. The 
key ingredient that allows dark stars to grow so much larger than ordinary 
fusion powered Population III stars is the fact that dark stars are so much 
cooler. Ordinary Pop 111 stars have much larger surface temperatures in 
excess of 50,000K. They produce ionizing photons that provide a variety of 
feedback mechanisms that cut off further accretion. McKee & Tan (2008) 
have estimated that the resultant Pop III stellar masses are ~ 14OM0. Dark 
stars are very different from fusion-powered stars, and their cooler surface 
temperatures allow continued accretion of baryons all the way up to enormous 
stellar masses, > lO^M©. 

WIMP annihilation produces energy at a rate per unit volume 



where is the WIMP number density, is the WIMP mass, and is 
the WIMP energy density. The annihilation products typically are electrons, 
photons, and neutrinos. The neutrinos escape the star, while the other anni- 
hilation products are trapped in the dark star, thermalize with the star, and 
heat it up. The luminosity from the DM heating is 



where fg is the fraction of the annihilation energy deposited in the star (not 
lost to neutrinos) and dV is the volume element. We take /q = 2/3 as is 
typical for WIMPs. 

Typically (100 — 10^^)71^0 of dark matter (up to ~ 1% of the halo mass) 
must be consumed by the star in order for large SMDS masses ~ lO^M© 
to be reached. We will consider two different scenarios for supplying this 
amount of dark matter: 



Qdm = nl{(Tv)m 



(2) 




(3) 
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1) Extended Adiabatic Contraction, labeled "without capture" below. In this 
case, dark matter is supplied by the gravitational attraction of the baryons in 
the star. The amount of DM available for DM annihilation due to adiabatic 
contraction may be larger than our previous estimates which were based on 
the assumption that DM halos are spherical. In a non-spherical DM halo, 
the supply of DM available to the star can be considerably enhanced, as we 
discuss in more detail in §2.21 In any case, this mechanism relies solely on 
the particle physics of WIMP annihilation and does not include capture of 
DM by baryons (discussed below). 

2) Extended Capture, labeled "with capture" below. Here, the star is 
initially powered by the DM from adiabatic contraction (AC), but the AC 
phase is taken to be short ~ 300, 000 years; once this DM runs out the star 
shrinks, its density increases, and subsequently the DM is replenished inside 
the star by capture of DM from the surroundings (Freese et al. 2008b; locci 
2008) as it scatters elastically off of nuclei in the star. In this case, the 
additional particle physics ingredient of WIMP scattering is required. This 
elastic scattering is the same mechanism that direct detection experiments 
(e.g. CDMS, XENON, LUX, DAMA) are using in their hunt for WIMPs. 
In previous work (Freese et al. 2008a; Spolyar et al. 2009), we assumed 
minimal capture, where DM heating and fusion contributed equally to the 
luminosity once the star reached the main sequence. Here, we consider the 
more sensible case where DM heating dominates completely due to larger 
ambient DM density, and the star can grow to become supermassive. 

Supermassive dark stars can result from either of these mechanisms for 
DM refueling inside the star. The SMDS can live for a very long time, 
tens to hundreds of million years, or possibly longer (even to today). We 
find that ~ lO^M© SDMSs are very bright ~ 3 x IO^Lq which makes them 
potentially observable by James Webb Space Telescope (JWST). We also 
note that SMDS may become even more massive a) if they form in larger 
haloes or b) the DM haloes in which they initially form merge with other 
haloes so that there is even more matter to accrete (n.b. alternatively these 
mergers may remove the DS from their high DM homes and stop the DM 
heating). For example, if dark stars form in lO^M© haloes, then they could 
in principle grow to contain all the baryons in the halo, i.e. M^, > 10 ''M©. 
Since the luminosity scales cis oc these SMDS would be even brighter, 

> IO^^Lq and are hence even better candidates for discovery in JWST. 
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Once the SMDS run out of DM fuel, they contract and heat up. The core 
reaches lO^K and fusion begins. As fusion-powered stars they don't last very 
long before collapsing to black holes. Again, this prediction is different from 
standard Pop III stars, many of which explode as pair-instabilty supernovae 
(Heger & Woosley 2002) with predicted even/odd element abundance ratios 
that are not (yet) observed in nature. The massive black holes (BH) remnants 
of the SMDSs are good candidates for explaining the existence of 10^ Mq 
BHs which are the central engines of the most distant {z > 5.6) quasars in 
the Sloan Digital Sky Survey (SDSS) (Fan et al. 2001, 2004, 2006) |. 

The idea of supermassive DS and the resultant > lO^M© BH was origi- 
nally proposed by Spolyar et al. (2009). Subsequently Umeda et al. (2009) 
took their existing stellar codes and added DM annihilation to allow the mass 
to grow. They started from Population III stars in which fusion was already 
present, assuming they then encounter a reservoir of DM. Our work differs 
in that we start from the very beginning with collapsing protostellar clouds 
that transition into dark stars, which can be DM powered for millions to 
billions of years before fusion ever begins. Our SMDSs are truly primordial 
supermassive stars. 

Begelman (2009) presents another alternative for the formation of su- 
permassive stars: rapid accretion onto stars which already have hydrogen 
burning in them. His "quasistars" , another possible route to large BH, are 
quite different from the SDMS discussed in this paper. 

Various other authors that have explored the repercussions of DM heating 
in the first stars (Taoso et al. 2008; Yoon et al. 2008; Ripamonti et al. 2009; 
locco et al. 2008). The possibility that DM annihilation might have effects 
on today's stars was initially considered in the '80s and early '90s (Krauss et 
al. 1985; Bouquet & Salati 1989; Salati & Silk 1989; Graff & Freese 1996) and 
has recently been studied in interesting papers by Moskalenko & Wai (2007), 
Scott et al. (2007), Bertone & Fairbairn (2007) and Scott et al. (2009). 
Other constraints on DS will arise from cosmological considerations. A first 
study of their effects (and those of the resultant MS stars) on reionization has 
been done by Schleicher et al.(2008, 2009) and further work in this direction 
is warranted. 

In this paper we examine the SMDS that result from the two mechanisms 
discussed above for DM refueling inside the star. In Section II we discuss the 
procedure for calculation of models; in Section HI we present results; and in 

^We thank N. Yoshida for pointing this out to us 
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Section IV we end with a discussion. 



2 Structure and Evolution of the Dark Star 

DM heating is very different from fusion. In order to overcome the Coulomb 
barriers between nuclei, fusion requires very high temperatures and densities 
in the star. Fusion is not very efficient in that only < 1% of the nuclear mass 
is converted to heat. WIMP annihilation, on the other hand, takes place at 
high dark matter densities regardless of the temperature. It is almost 100% 
efficient since 0(1) of the WIMP mass is converted to useful energy. Thus, 
in the evolution of the first protostars, DM heating becomes important early. 
Here we start the calculation when the DS is massive enough (3 Mq) so 
that it is in hydrostatic equilibrium and most of the hydrogen and helium is 
ionized. The contribution to DM luminosity is roughly constant as a function 
of radius throughout the DS, unlike fusion which takes place only at the (high 
temperature) core of the star. 



2.1 Basic Equations 

We use the numerical code previously discussed in detail by Freese et al. 
(2008a) and Spolyar et al. (2009). We make the assumption that a DS can 
be described as a polytrope 

P = (4) 

in hydrostatic equilibrium. Here P is the pressure, p is the density, and K 
is a constant. We solve the equations of stellar structure with polytropic 
index n initially 1.5, as appropriate for convective stars, and made a gradual 
transition to n = 3 as the star becomes radiative in the later phases. We 
require that at each time-step during the accretion process the star is in 
hydrostatic equilibrium, 

dP _ GMr 
dr 



-Pir)^ (5) 



where = 47rr^p(r), and Mr is the mass enclosed in a spherical shell of 
radius r. The equation of state includes radiation pressure, 

WMIW^i 4, (6) 
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where is the Boltzmann constant, m„ is the atomic mass unit and /i = 
0.588 is the mean atomic weight. The opacity is obtained from a zero metal- 
hcity table from OPAL (Iglesias & Rogers 1996) supplemented at low tem- 
peratures by opacities from Lenzuni et al. (1991) for T < 6000K. The further 
assumption is made that the radiated luminosity of the star L^, is balanced 
by the rate of energy output by all internal sources, Ltot, as described below 
in 3231 

= AnRlaBT^ff = Ltot- (7) 

where Tg// is the surface temperature, and i?* is the total radius. 

Starting with a mass M and an estimate for the outer radius i?*, the 
code integrates the structure equations outward from the center. The total 
rate of energy production Ltot is compared to the stellar radiated luminosity, 
as in equation ([7]) and the radius is adjusted until the condition of thermal 
equilibrium is met (a convergence of 1 in 10^ is reached). 

2.2 Dark Matter Densities 

We now describe the two different mechanisms for supplying the DM density. 

Extended Adiabatic Contraction: As the baryons start to collapse into a 
protostellar cloud at the center of the DM halo, the DM responds to the 
changing gravitational potential well and falls in as well. As described in our 
previous work (Spolyar et al. 2008), we will use adiabatic contraction (AC) 
to describe this increase in DM density. For the case of spherical haloes, we 
previously found, by performing exact calculations for comparison (Freese 
et al. 2009), that the simple Blumenthal method (Blumenthal et al. 1985; 
Barnes & White 1984; Ryden & Gunn 1987) gives reliable results for the 
final DM densities up to an unimportant factor of two; others confirmed this 
conclusion (Natarajan et al. 2009; locco et al. 2008; Sivertsson & Gondolo 
2010). Using this simple approach during the AC phase, we found that 

~ 5GeV/cm^(?T,/i/cm^)°'^-'^ where Uh is the gas density. These are the 
values we will use during AC. 

In our previous work, we probably underestimated the lifetime of the 
DM inside the star due to AC. In our previous work we treated the DM halo 
as spherical and ran up the DS mass to the point where the DM initially 
inside the star was entirely consumed by annihilation. The DS mass at this 
point is O(1O^)M0 after a lifetime of ~ 300, 000 years, and the amount of DM 
consumed has only amounted to ~ IMq. In a spherical DM halo, the orbits of 
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DM particles are planar rosettes (Binney & Tremaine 2008) conserving energy 
as well as all three components of angular momentum; consequently the 
central hole (or "empty loss cone") that results from DM annihilation cannot 
be repopulated once it is depleted. (Note that although DM annihilation 
creates a central hole in the DM density, the entire region is filled with 
baryons and hence the potential is stable.) However, it is well known that 
DM halos formed in hierarchical structure formation simulations are not 
spherical but are prolate-triaxial (Bardeen et al. 1986; Barnes & Efsthatiou 
1987; Frenk et al. 1988; Dubinski & Carlberg 1991; Jing & Suto 2002; Bailin 
& Steinmetz 2005; AUgood et al. 2006) with typical axis ratios of (short- 
axis)/(long-axis) ~ 0.6 — 0.8. In triaxial potentials, the orbits do not conserve 
angular momentum. In particular there are two families of "centrophilic 
orbits" (box orbits and chaotic orbits) which oscillate back and forth through 
the potential and can travel arbitrarily close to the center (Schwarzschild 
1979, Goodman & Schwarzschild 1981; Gerhard & Binney 1985; de Zeeuw 
1985; Schwarzschild 1993; Merritt & Pridman 1986; Merritt & Valluri 1996). 
Unlike an orbit in a spherical potential which has a constant pericenter radius 
(the distance of closest approach to the center of the potential), the pericenter 
radius of a centrophilic orbit varies over time extending from rperi = [0, rmax] 
where rmax is sometimes referred to as the "throat" of the orbit (Gerhard 
& Binney 1985). As any one particle traverses the center of the dark star, 
it may indeed be removed from the pool by annihilation. However, it was 
unlikely to pass through the star (of radius r*) on its next orbit anyhow since 
in general Vmax >> t*. Instead, a particle on a different "centrophilic orbit" 
enters the star for the first time maintaining the steady state central DM 
density cusp. Unlike in the case of a spherical DM halo, where annihilation 
steadily depletes the central density cusp, in a triaxial halo there is a high 
probability for a particle on a centrophilic orbit to pass through the center 
for the first time on any given orbital crossing. This is particularly true in 
potentials with central point masses which are dominated by chaotic orbits 
and are therefore ergodic (Merritt & Valluri 1996). Hence, the central DM 
density can remain much higher than we previously expected. 

The dynamics of the refilling of the central "loss cone" in the case of 
spherical and non-spherical coUisionless systems has been studied previously 
in the context of capture of stars by a central black hole (Gerhard & Binney 
1985; Magorrian & Tremaine 1999; Merritt & Poon 2004). The details of 
filling rate for the specific case of the DS will be presented separately in a later 
paper; however these previous papers showed that in non spherical systems 
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the loss cone could remain full for a factor of 2 to 10^ times longer than 
in the spherical case depending on whether the potential was axisymmetric 
(Magorrian & Tremaine 1999) or triaxial (Merritt & Poon 2004), respectively. 
Since DM halos are known to be triaxial this suggests that the duration for 
which the central hole remains full (i.e. has orbits with the low angular 
momentum necessary for annihilation) can increase from 300,000 years to as 
much as 3 X 10^ yrs possibly allowing the DS to be detectable by JWST. 

These more optimistic estimates require that a significant fraction of the 
orbits in these early DM halos are chaotic and boxlike. One important poten- 
tial concern with assuming conditions in a triaxial halo is that several studies 
have shown that the growth of central baryonic components tend to make 
DM halos more axisymmetric than in purely dissipationless simulations (Du- 
binski 1994, Evrard et al. 1994, Merritt & Quinlan 1998, Kazantzidis et al. 
2004, Debattista et al. 2008, Tissera et al. 2009), and axisymmetric models 
are generally not expected to contain centrophilic orbit families. However, 
Valluri et al. (2010) recently showed that when a compact central baryonic 
component is grown adiabatically inside a triaxial DM halo, the final halo 
that results from the adiabatic growth of such a baryonic component looks 
nearly oblate, yet its orbit population can contain a significant fraction of 
centrophilic orbits, since box orbits preferentially deform their shapes rather 
than converting to centrophobic tube families. Furthermore they showed that 
a significant fraction of the orbits (both box orbits and tube orbits which were 
traditionally thought to avoid the center) experience strong chaotic scatter- 
ing, a mechanism that could drive them close to the center. They found that 
for a fixed ratio of the mass of the baryonic component to mass of the DM 
halo (M*/Mdm = 10~^), the smaller the radius of the baryonic component 
(r^), relative to the virial radius of the DM halo (ryir) the greater was the 
fraction of chaotic orbits. The most compact baryonic component studied 
by them had r^/r^u = 4.6 x 10^'^, which is significantly larger than that for 
our fiducial SMDS for which r*/rvir = 2 x 10~^. The compactness of the 
baryonic component in the DS relative to its halo is important since Valluri 
et al. (2010) showed that when the central component became very com- 
pact, orbits that were previously thought to be immune to becoming chaotic 
(the "long-axis tubes" which are important in prolate DM halos) also be- 
come chaotic. Valluri & Preese (2010) are currently computing the rate at 
which the "loss cone" at the center of the SDMS will be refilled for a range 
of possible halo and DS masses. 

It is interesting to speculate that the Initial Mass Function of the first 
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stars may be determined by the cutoff of the DM supply, which will take 
place at different SMDS masses in different haloes, depending on the details 
of the cosmological merger history. As the SMDS mass becomes a significant 
fraction of the halo mass (~ 5 — 10%), it can significantly affect the shape of 
the halo, causing it to become more axisymmetric (Merritt & Quinlan 1998) 
and thereby potentially cutting off the DM supply; on the other hand, in 
the meantime the halo is growing larger due to mergers, which will replenish 
the population of radial orbits. Numerical simulations with better resolution 
than currently possible will be necessary to address these questions. 

In the meantime, in our case without capture, we assume that the required 
DM is present and allow the stellar mass to grow to the point where most of 
the baryonic content of the initial halo is inside a single SMDS. 

The amount of dark matter required inside the star to sustain long enough 
DM heating to reach a lO^M© SMDS is stiU smaU, ~ lOOM© for accretion 
rate M = lO^^Mo/yr and ~ lO^M© for accretion rate M = IQ-^Mq/jt, 
out of a total IO^Mq halo. In the code we accomplish this by not removing 
annihilated DM from the pool. More precise studies must be performed later 
in which we follow individual particle orbits in triaxial potentials to better 
determine the precise DM density at any one time. 

Extended Capture: In our model labeled "with capture", we assume (as 
in our previous papers) that the DM due to adiabatic contraction runs out 
in ~ 300, 000 yrs. For a while DM heating becomes unimportant and the 
DS has to contract to maintain pressure support. Then DS is refueled in the 
later stages due to capture of further DM from the ambient medium. This 
refueling requires an additional piece of particle physics: scattering of DM 
off the nuclei inside the star. 

Some of the WIMPs from the ambient medium that have orbits passing 
through the star will eventually be captured and sink to the center, where 
they can annihilate efficiently The capture process is irrelevant during the 
early evolutionary stage of the DS, since the baryon density is not high 
enough at that point, leading to very small scattering probabilities. However, 
once the DS approaches the main sequence, the baryon densities become 
high enough for substantial capture to be possible. This mechanism was first 
noticed by Freese et al. (2008b) and locco (2008). 

In our previous work (Spolyar et al. 2009) we investigated a 'minimal 
capture' case which did not cause the DS to grow much more massive than 
the original case without capture ~ lOOOM©; but we stated our intention to 
work out the case of a more substantial background DM density in which 
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case the DS would end up supermassive. This is what we investigate here. 
The capture rate is sensitive to the product of two uncertain quantities: 
the scattering cross section of WIMP interactions with the nuclei cxc and 
the ambient DM density p^. Since the capture mechanism depends only 
on the product of these two quantities, one can interchangeably vary either 
of these. For illustration purposes we will fix (Xc = 10~^^ cm^ and vary 
= (10^° — 10^'*)GeV/cm^. The latter quantity is the largest reasonable 
amount based on our results for AC at the DS surface; and the former is the 
"minimal capture" value considered in all our previous papers. Our fiducial 
cross section is just below the experimental bound for spin- dependent (SD) 
scattering; the bound on spin-independent (SI) scattering is much tighter: 
ac,si < 3.8 X 10"^^ cm2 for = lOOGeV (Ahmed et al. 2009). We will 
show that capture can produce sufficient DM in the star to keep DM heating 
alive for a long time. The details of our procedure for including capture have 
previously been presented in Spolyar et al. (2009) and will not be repeated 
here. 

2.3 Energy Sources: 

There are four possible contributions to the DS luminosity: 

grav ~ ^nuc \ -'-'cap (8) 

from DM annihilation, gravitational contraction, nuclear fusion, and cap- 
tured DM respectively. The heating due to DM annihilation in Eqns. ([2]) 
and (E]) dominates from the time of DS formation until the adiabatically 
contracted DM runs out. As described previously, in our "without capture" 
models this stage never ends due to extended adiabatic contraction. In the 
models "with capture", on the other hand, we take this phase to end after 
~ 300,000 years, so that the DS has to contract in order to maintain pres- 
sure support. The contribution Lgrav due to gravitational energy release is 
calculated as in Spolyar et al. (2009). As the DS contracts, the density 
and temperature increase to the point where nuclear fusion begins. We in- 
clude deuterium burning starting at T ~ lO^K, hydrogen burning via the 
equilibrium proton-proton cycle (Bahcall 1989), and helium burning via the 
triple-alpha reaction (Kippenhahn & Weigert 1990). During the later stages 
of the pre main sequence evolution in the cases "with capture", the DS be- 
comes dense enough to capture DM from the ambient medium via elastic 
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scattering. Already before fusion can begin, and possibly again after the 
onset of fusion, captured DM can provide an important energy source with 
accompanying luminosity 

Leap = '^m^T cap = '2rnjQ J dV pl^p{av) / (9) 

and again /q = 2/3. 



3 Results of Stellar Structure Analysis 

Using our polytropic model for dark stars, we have started with 3Mq stars 
and allowed baryonic matter to accrete onto them until they become super- 
massive with > IO^Mq. We display results for the case without capture 
(but with extended adiabatic contraction) as well as the case with capture for 
a variety of WIMP masses = 10 GeV, 100 GeV, and 1 TeV. We have run 
models for a variety of accretion rates of baryons onto the star including con- 
stant accretion rates of M = 10^"'^, 10^^, 10~^MQ/yr. We will present results 
for M = lO~^M0/yr, which is approximately the average rate calculated by 
Tan & McKee (2004) and by O'Shea & Norman (2007). 

Our stellar evolution results can be seen in the Hertzsprung-Russell di- 
agram of Figure 1 for the case of a IO^Mq halo. The dark star travels up 
to increasingly higher luminosities as it becomes more massive due to accre- 
tion. We have labeled a sequence of ever larger masses until all the baryons 
(15O,OOOM0) in the original halo are consumed by the SMDS. As the mass 
increases, so do the luminosity and the surface temperature. In the cases 
"without capture," the radius increases continuously until all the baryons 
have been eaten. In the cases "with capture" , we have taken the (overly con- 
servative) assumption that the DM from adiabatic contraction is depleted 
after ~ 300, 000 yrs as in our earlier papers; then the luminosity plateaus for 
a while while the DS contracts until eventually it is dense enough to capture 
further DM. 

We note that, for the case "without capture", the H-R diagram is un- 
changed by varying the accretion rate: only the time it takes to get from one 
mass stage to the next changes, but the curves we have plotted apply equally 
to all accretion rates. Similarly, given m^^, the following quantities are the 
same regardless of accretion rate: i?*,Te//,pc, and Tc- 
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In a beautiful paper, Hoyle & Fowler (1963) studied supermassive stars 
in excess of and found results germane to our work. They treated 

these as n = 3 polytropes (just as we do) dominated by radiation pres- 
sure, and found the following results: ~ W^{M^/MQy/^{Tc/10^K)-^cm, 
U/Lq ~ lO^M./Mo, and T^jf ~ 10%Tc/10^Ky/'^K. While some of the 
details of their calculations differ from ours, taking the central temperature 
appropriate to DS in the above relations roughly reproduces our results (to 
0(1)). For example, by using the temperature appropriate to dark stars with 
extended AC (Tc ~ lO^K) rather than the much higher central temperature 
(Tc > lO^K) appropriate to nuclear power generation, the above relations 
show that DS have much larger radii and smaller surface temperatures than 
fusion powered stars. We wish to draw particular attention to the fact that 
luminosity scales linearly with stellar mass, and is independent of power 
source. 

Figure 2 plots the H-R diagram "with capture" for a single WIMP mass 
of 100 GeV, for M = lO'^ Mq/jt, and for = lO'^^cm^, but for a variety 
of ambient densities ranging from = (10^° — 10^^)GeV cm~^. The latter is 
the density one finds due to adiabatic contraction at the photosphere of the 
DS, and seems the largest sensible starting point for the value of the ambient 
density. Our previous paper (Spolyar et al. 2009) considered the case of 
minimal capture with 10^°GeV/cm^, which was artificially chosen so that the 
growth of the DS ceases at 68OM0, the radius shrinks, and then fusion and 
DM heating play equal roles. For ambient densities below 5 x 10^*^GeV/cm^, 
the DS mass growth shuts off well before the star becomes supermassive for 
the following reason. The cases with capture all take place at higher stellar 
densities than the cases without; the density must be high enough to be 
able to capture WIMPs. Consequently the surface temperature is larger and 
accretion is shut off more easily by radiation coming from the star. The 
case of ambient density 10^°GeV/cm^ (from our previous paper) is a very 
carefully chosen (dehcate) situation. On the low side of this density, DM 
heating is completely irrelevant and fusion tells the whole story; on the other 
hand, for any density > few x 10^'' GeV/cm^, DM heating is so dominant 
over fusion that the DS can just continue growing in mass. At these higher 
densities the surface temperature never becomes hot enough (pa 100, 000 K) 
for feedback effects from the star to cut off accretion. Between 50,000 K and 
100,000 K feedback effects are included, and they act to reduce the accretion 
rate, but they never shut it off entirely for densities above 5 x 10^°GeV/cm^. 
In reality a star that is moving around can sometimes hit pockets of high p^ 
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(where it is DM powered and grows in mass) and sometimes hit pockets of 
low (where fusion takes over). As long as the ambient density remains at 
least this large, the star can reach arbitrarily large masses and eat the entire 
baryonic content of the halo. 

As described previously, the capture mechanism depends on the product 
of scattering cross section times ambient density, cycPxi rather than on either 
of these quantities separately. Hence our current discussion could trade off 
ambient density vs. cross section. For example, the product is the same for 
p-^ = 5 X 10^° GeV/cm^ and cTc = 10~^^ cm^ as it is for p^^ = lO^^GeV/cm^ 
and (Tc = 5 X 10~^^ cm^. Thus, for the highest reasonable ambient density, 
the scattering cross section can be several orders of magnitude lower than 
the experimental upper bound for spin- dependent (SD) scattering and still 
provide substantial capture in DS. While the required cxc is ruled out (Ahmed 
et al. 2009) for SI scattering at = lOOGeV, it is below the bounds at low 
masses ~ 10 GeV and in this case can lead to significant DM capture in 
the stars. 

Above ~ IOOMq, one can see that the stellar luminosity scales 
and is the same for all models for a given stellar mass; this statement is 
essentially true for all stars no matter the power source. The reason is that at 
these masses, the star is essentially radiation pressure supported throug houli. 
This same scaling in supermassive stars was already noticed by Hoyle and 
Fowler (1963). The luminosity essentially tracks (just below) the Eddington 
luminosity which scales as L cx M^,. 

The curves with higher values of WIMP mass lie to the left of the 
curves with lower m^. This can be understood as follows. The DM heating 
rate in Eqn.(l) scales as Q cc p'^/m^. Hence to reach the same amount of 
heating and achieve the same luminosity, at higher the DS must be at 
higher WIMP density, i.e., the stellar radius must be smaller, the DS is hotter, 
and the corresponding surface temperature Te// is higher. Also, for higher 

the amount of DM in the star is smaller since the star is more compact for 
the same number of baryons but oc n^'^ where n is the hydrogen density. 

Tables 1 and 2 show various stellar properties for a DS that forms in a 
IO^Mq dm halo, as the star evolves to higher mass for the case of = 100 
GeV, for the two cases "without" and "with" capture respectively. While 
the DM density is a gently decreasing function of radius for the case without 

■^There is a slight deviation for the case without capture, where the star is stih 

only 78% radiation pressure supported with the remaining pressure due to gas. 
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capture, it is extremely sharply peaked at the center of the DS for the case 
"with capture". 

One can see that, in the case "without capture", ~ 10'* Mq of dark matter 
must be annihilated away in order for the DS to reach lO^M© for accretion 
rate M — 10~^Mq. The time to reach this mass is ~ lOOMyr. For an 
alternative faster accretion rate M — 10~^Mq, a smaller amount of DM 
must be annihilated away, ~ IOOMq; then it takes the DS ~ 10 Myr to reach 
~ IO^Mq. The caveat is that the DM orbits must continue to penetrate into 
the middle of the DS for this length of time in order for the DM abundance 
and consequent heating to continue; it is the DM heat source that keeps the 
DS cool enough to allow it to continue to accrete baryons. Additionally, the 
assumption that baryons continue to accrete onto the DS must continue to 
hold. Yet, in the time frame required, the original IO^Mq halo will merge with 
its neighbors, so that both the baryon and DM densities are disturbed. These 
mergers could effect the DS in one of two ways: either they provide more 
baryons and DM to feed the SMDS so that it ends up being even larger, or 
they disrupt the pleasant high DM environment of the SMDS so that it loses 
its fuel and converts to an entirely fusion powered star. Continued growth 
of the DS is quite plausible since simulations with massive BHs in mergers 
show that they prefer to sit close to the center of the density distribution or 
find the center in a short time after the merger. 

Someday detailed cosmological simulations will be required to answer 
this question. Individual DS in different haloes may end up with a variety 
of different masses depending on the details of the evolution of the haloes 
they live in. The case studied in this paper is clearly a simplistic version of 
the more complicated reality, but illustrates the basic idea that supermassive 
stars may be created by accretion onto DS, either with or without capture. 

Tables 1 & 2 present models for DSs which form in DM halos of mass 
for a 10^ Mq, but SMDS could also form in a variety of halo masses with 
different final stellar masses. For example, a hydrogen/helium molecular 
cloud may start to contract in a lO^M© halo and produce a DS. Here the 
situation is more complicated. The virial temperature of the halo exceeds 
lO'^K, the surface temperature we have found for a DS in equilibrium. Hence 
it is not clear how accretion onto the DS will proceed. This is the subject 
of future work. The accretion is expected to be faster due to the increased 
ambient temperature. We extended our models to IO^Mq halos in which 
a more extended period of accretion can lead to SMDS with even larger 
masses. Tables 3 & 4 show examples of "without" and "with" capture cases 
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Table 1: Properties and evolution of dark stars for = 100 GeV, 
M — 10~^ Mq/jt for the case without capture but with extended adiabatic 
contraction. The DM halo was taken to be at a redshift of 20 with a concen- 
tration parameter of 3.5 and with a mass of IO^Mq. Shown are the stellar 
mass M*, the DS luminosity L^,. the stellar radius R^, the surface temper- 
ature Tgff, the central baryon density Pc, the central temperature T^, the 
amount of DM in the DS (due to both adiabatic contraction and cap- 
ture) , the central WIMP density p^_c, and the amount of DM consumed by 
the DS Mauu- 
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0.9 
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21 


5.3 


8.4 X 10-^ 


86 
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26 


16.5 


4.1 X 10"^ 


46 


31.2 


1.6 X 10~« 


10750 
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Table 2: Properties and evolution of dark stars for case "with capture", for 
— 100 GeV, M — lO'^Mg/yr, and product of scattering cross section 
times ambient DM density acPx = 10~^^cm^ x 10^^ GeV/cm^. The Halo 
has the same parameters as in Table 1. The quantities tabulated are the 
same as in Table 1. The double horizontal line delineates the transition from 
adiabatically contracted DM to captured DM once the DS reaches ~ IOOOMq 
(after this point, the DM from AC has been annihilated away). 
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respectively for IO^Mq halos. The baryonic mass in the halo is 1.5 x 10"Mo. 
Potentially all of this mass could go into the SMDS. Then once the DM runs 
out it becomes an enormous Pop 111 fusion powered star, which soon burns 
out and becomes a black hole with mass > lO^Mo. Such a large BH at early 
times would clearly help to explain the many large BH found in the Universe 
at early times and today. In addition, since oc M*, we can predict that the 
luminosity of the ~ lO'^ Mq SMDSs would be ~ IO^^Lq, even easier 
to detect than a 10^ SMDS. 

General Relativistic Instability: The pulsational stability of supermassive 
stars is an interesting issue. They are radiation-pressure dominated with 
adiabatic index close to 7 = 1 + 1/n = 4/3, the value that yields neutral 
stability to radial pulsations for Newtonian bodies with no rotation. Indeed 
general relativistic corrections (which scale as GM^/R^) act in the direction 
of destabilizing stars and are particularly important for high mass stars. 
Fowler (1966) examined the stability of supermassive stars using polytropes 
with n — ?> (see Wagoner 1969 for a review). Fowler found that, for the case 
of no rotation, radial oscillations become dynamically unstable and prevent 
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Table 3: Properties and Evolution of dark stars for rriy^ = 100 GeV, 
M — 10"^ Mq/jic for the case without capture but with extended adiabatic 
contraction. The DM halo was taken to be at a rcdshift of 15 with a concen- 
tration parameter of 3.5 and with a mass of lO^M©. The quantities tabulated 
are the same as in Table 1. 
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standard stars more massive than lO^M© from reaching a phase of hydrogen 
burning before collapse. Yet he also found that a small amount of rotation 
can stabilize the stars, so that rotating stars as heavy as 10® Mq could be 
stable en route to reaching hydrogen burning. 

In the case of DS, stabihty to radial pulsations is much easier to achieve. 
DS have much larger radii and lower temperatures than fusion powered stars, 
so that the GR corrections ~ GM^/R^ are much smaller. The upper limit on 
the allowed stellar mass will be larger. In any case, SMDSs are undoubtedly 
rotating, so that very large stable masses can be achieved (even in the case of 
rotating ordinary stars, the mass limit is 10® M©). In the future we suggest 
a stability analysis of our models. 

4 Detectability with JWST 

Wc discuss the capabilities of JWST to discover dark stars, following the 
properties of the telescope described by Gardner et al. (2006, 2009). The 
telescope is designed to be diffraction limited at a wavelength Ao^^ — 2/xm. 
The Near Infrared Camera (NIRCam) will operate in the wavelength range 
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Table 4: Properties and evolution of dark stars for case "with capture", for 
— 100 GeV, M — 10"^ Mq/yt, and product of scattering cross section 
times ambient DM density acPx = lO'^^cm^ x lO^'^ GcV/cm\ The DM 
halo has the same parameters as in Table 3. The quantities tabulated are 
the same as in Table 1. The double horizontal line delineates the transition 
from adiabatically contracted DM to captured DM once the DS reaches ~ 
4 X 10* (after this point, the DM from AC has been annihilated away). 
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A = (0.6 — 5)/im and the Mid- Infrared Camera (MIRI) will operate in 
the wavelength range A = 5 — 27/im. In an exposure of duration lO'^s, 
NIRCam will have a limiting sensitivity of 11.4 nJy (InJy = 1 x 10~^^ 
ergs s~^ cm~^ B.z~^) in the 2/xm band, and 13.8nJy in the 3.5/im band and 
MIRI will have a sensitivity of 700/iJy in the A = lOfim band (in all cases 
limiting sensitivities are for a S/N=10). With longer exposure times the 
limiting flux detectable will scale as \/t for the same S/N. DS will be char- 
acterized by black body spectra with surface temperatures Te// ^ 5 x lO^K. 
In addition, DS are also predicted to have hydrogen lines. 

We determine the detectability of dark stars located at various redshifts 
z= 5, 10, and 15 using the standard Planck spectrum of blackbody with 
surface temperature T^ff and radius i?* (for DSs similar to those from Tables 
1-4) in a cosmology with Hq = 74, I^a = 0.71, ^Im = 0.29. 

Figure |3] shows the observed black body flux distribution of two SMDSs 
formed at z = 15 for a WIMP mass = lOOGeV for the case of extended 
AC (without capture). The star in the left panel is formed in a IO^Mq 
halo and the star in the right panel is formed in a 10^ Mq halo and their 
stellar (baryonic) masses are 1.7 x 10^ Mq and 1.5 x 10^ Mq respectively. 
Curves are shown assuming the SMDS formed at ^ = 15 and survived to 
various redshifts, at which it is still producing blackbody radiation. The 
1.7 X 10^ M0 star (left panel) will be detectable by JWST (NIRCam) in an 
exposure of a million seconds, but only if it survives intact till 2; = 10. The 
1.5 X 10'' Mq star (right panel) will be detectable even in a shorter lO^s 
exposure even at z = 15 in both the 2/i and 3.5/i bands. The star on the 
right may be marginally detectable in a million second exposure in the 10// 
band of MIRI. The relative flux levels in the three different bands will be 
important for distinguishing these objects from galaxies. 

The curves are not corrected for Ly-a absorption by the IGM but the 
red vertical lines show the location of the 12 16 A line redshifted from the 
rest-frame wavelength of the star at each of the three redshifts. Flux at 
wavelengths to the left of the redline at each redshift is expected to be ab- 
sorbed to some extent by the IGM. Since the surface temperatures of our 
stars are ~ lO^K, the majority of the Lya absorption {Xrest = 1216A) is 
expected to occur at wavelengths shorter than that at which the peak flux 
is emitted. We note that for a DS at z = 15 the Ly a absorption line lies at 
1.94/i - roughly in the middle of the NIRCam 2fi band. In this case the flux 
in this band will be reduced by about a factor of two but will still be well 
above the detection limit (in lO^s). A detailed calculation of the absorption 
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by the intergalactic medium is outside the scope of the current paper. 

We studied numerous other cases without capture as welh a variety of 
WIMP masses (lOGev -2TeV) as well as various formation redshifts (zform = 
10 — 20). We found that DSs with masses up to M* ~ lO^M© forming in 
10^ Mq halos at 2 = 20 and shining at that redshift will in general not 
be detectable for any value of (they will become detectable only if they 
survive to and shine at much lower redshifts). 

The smaller DS with = 8OOM0 discussed in our previous work (Freese 
et al. 2008a; Spolyar et al. 2009) will not be visible in JWST (see also Scott 
et al 2010); they are several orders of magnitude below the detection limit of 
JWST (in lO'^s). 

Figure m shows the observed black body flux distribution of two dark stars 
formed at z = 15 in halos of two different masses for the case "with capture" . 
Since DS formed via capture are smaller (in radius) and hotter (than DS 
formed via extended AC without capture), their peak wavelength tends to 
shift to lower wavelengths, in some cases out of the range detectable by 
JWST. lO^M© stars in IO^Mq halos are only detectable if they survive until 
z = 5 at which time they could be detectable in a long (lO^s exposure). On 
the other hand 10'^ stars formed in lO^M© halos will be easily detectable 
even in an exposure of lO^s all the way out to z = 15. For most other WIMP 
masses and formation redshifts DSs formed in lO^M© halos via capture are 
below the detection hmit of JWST. 

The prospect of detecting SMDS in JWST and confirming the existence 
of a new phase of stellar evolution is exciting. In the most optimistic cases, 
detection in a number of different wavelength bands could be used to obtain 
a spectrum and differentiate these dark stars from galaxies or other sources. 

5 Concluding Remarks 

Using our polytropic model for dark stars, we have considered accretion of 
baryonic matter onto the DS as they become supermassive, > IO^Mq. 
Such large masses are possible because the dark star is cool enough (as long 
as it is powered by DM) so that radiative feedback effects from the star do 
not shut off the accretion of baryons, as long as it is powered by DM. We 
considered two different scenarios for supplying the required amount of dark 
matter: 

1) The Case of Extended Adiabatic Contraction, labeled "without 



21 



capture" in the figures. In this case dark matter is supplied by the 
gravitational attraction of the baryons in the star. In triaxial haloes DM 
orbits are quite complex and the DM in the core is harder to deplete than 
previously estimated. This case does no^ include any captured DM, and relies 
solely on the particle physics of WIMP annihilation. To grow to a lO^M© 
SMDS in a IO^Mq halo, or to grow to a WMq SMDS in a lO^M© halo, the 
amount of DM consumed can be as much as ~ 1% of the total DM in the 
halo (depending on the accretion rate). This amount is not unreasonable, 
since Valluri etal (2010) found that the fraction of box and chaotic DM orbits 
is as high as 85% in triaxial haloes and remains over 10% when a significant 
compact baryonic component causes the halo to become axisymmetric at 
small radii. Future work will be required to accurately obtain WIMP orbits, 
densities, and timescales (work in progress). For now we took the simplistic 
approach of using our previous prescription for adiabatic contraction in a 
spherical potential but not removing the annihilated DM. 
2) The Case of Extended Capture, labeled "with capture" in the 
figures. Here the original DM inside the star from adiabatic contraction is 
assumed to be depleted after ~ 300,000 yrs, then the star begins to shrink 
somewhat, and capture of DM from the surroundings takes place as it scatters 
elastically off of nuclei in the star. In this case the additional particle physics 
ingredient of WIMP scattering is required. 

In this paper we studied the formation of IO^Mq SMDS in lO^M© DM 
haloes and IO'^Mq SMDS in lO^M© haloes. These stars become very bright, 

~ (10^ — 1O^^)L0. Figure 1 shows the H-R diagram for a variety of WIMP 
masses, and follows the dark star as it chmbs up to ever higher masses. 
They live millions to billions of years, depending on the merger history with 
other haloes. Once the DM runs out, the SMDSs have brief lives as fusion 
powered Pop III stars before collapsing into > IO^Mq black holes, possible 
seeds for many of the big BH seen in the Universe today and at early times. 
A proper study of the final mass of the DS and resultant BH will depend 
on cosmological simulations. The original halo containing the DS will merge 
with other haloes. No one knows what exactly will happen to the DM density 
in the vicinity of the DS when this happens. The DS could end up even 
more massive. DS may also form in larger haloes that form at later times, as 
long as the baryonic content is still only H and He. Localized regions with 
this property could exist even at redshifts z < 7 (Furlanetto & Loeb 2005; 
Choudhury & Ferrara 2007). 

SMDS would make plausible precursors of the lO^M© black holes observed 
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at 2; > 6 (Fan et al. 2003) of intermediate mass black holes; of BH at the 
centers of galaxies; and of the BH inferred by extragalactic radio excess 
seen by the ARCADE experiment (SeifFert et al. 2009). In addition, the 
BH remnants from DS could play a role in high-redshift gamma ray bursts 
thought to take place due to accretion onto early black holes 

SMDSs could be detected by JWST for a variety of parameter ranges. 
These are extremely bright objects ~ (10^ — 10^^)Lq and yet are very 
cool T ~ 10, OOOK, so that their emitted light is in the wavebands detectable 
by JWST. The longer they live, the more easy they are to detect. Figures 
3 and 4 give examples of what one could look for in JWST. For the most 
optimistic cases, one could even test for the blackbody spectrum in a number 
of different wavebands. In principle hydrogen or helium lines could be found 
to complement the blackbody emission. If, in addition, someday high energy 
neutrinos are found to emanate from these stars, then it will be a clincher 
that DM annihilation took place inside the DS. 

It is interesting to speculate that the Initial Mass Function of Population 
III fusion powered stars may be determined by the cutoff of the DM supply, 
which may vary from one DS to another. Dark stars continue to accrete mass 
as long as the dark matter annihilation powers the star and keeps it cool 
enough. Once the DM fuel supply is exhausted, the star shrinks and heats 
up, fusion begins, and the mass growth of the star is quickly halted due to 
feedback from hot emitted photons. Hence the details of the cutoff of the DM 
supply may determine the sizes of Population HI stars entering the fusion 
era. The cutoff will take place at different DS masses in different haloes, 
depending on the details of the cosmological merger history. Different final 
DS masses may result for different individual DS depending on the evolution 
of their parent haloes. 
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Figure 1: Hertzsprung- Russell diagram for dark stars for accretion rate M = 
IQ-^Mo/yr. and a variety of WIMP masses as labeled for the two cases: 
(i) "without capture" but with extended adiabatic contraction (dotted lines) 
and (ii) "with capture" (solid lines). The case with capture is for product 
of scattering cross section times ambient WIMP density (TcP^ = 10~^^cm^ x 
lO^^GeV/cm^. Also labeled are stellar masses reached by the DS on its way 
to becoming supermassive. The final DS mass was taken to be 1.5 x 10^ Mq 
(the baryonic mass inside the initial halo), but could vary from halo to halo, 
depending on the specifics of the halo mergers. 
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Figure 2: Hertzsprung-Russell diagram for DS for the case "with capture" 
for 100 GeV WIMP mass and accretion rate M = lO~^M0/yr. The different 
curves are for a variety of ambient DM densities as labeled for scattering 
cross section = lO^^^cm^. The results depend only on the product cTcPx 
so the different curves could equivalently refer to different ctc for a given p^. 
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Figure 3: Black body spectra of two dark stars formed via extended adiabatic 
contraction ("without capture") for m^=100 GeV. Left panel: 1.7 x 10^ Mq 
SMDS in a 10*^ Mq halo. Right panel: 1.5 x 10^ M© SMDS in 10^ Mq halo. 
The black body flux is shown at z = 15 (formation redshift) and at 2 = 10 
and 5 (see line legends) assuming that the dark star survives till the lower 
redshifts. Blue dashes show sensitivity limit and bandwidth of NIRCam 2// 
(R=4) while the green dashes show the sensitivity limit and band width of the 
NIRCam 3.5/i (R=4) band. The upper and lower dashes show the sensitivity 
limits after exposure times of lO^s, lO^s respectively. The sensitivity of MIRI 
(10/i, R=5) is shown for exposure time of lO^s (orange dash). All sensitivities 
are computed assuming a S/N=10. The red vertical lines show the location 
of the 1216A line redshifted from the rest-frame wavelength of the star at 
each of the three redshifts. The observed flux to the left of the vertical lines 
will decrease relative to the black curves depending on the model assumed 
for IGM absorption up to the redshift of reionization. 
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Figure 4: Similar to Fig. |3]for dark stars formed "with capture". Left panel: 
1.7 X 10^ Mq SMDS formed in lO^M© halo {m^ = 50Gev). Right panel: 
1.7 X 10^ Mq SMDS formed in IO^Mq halo (m^ = lOOGev). 
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